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A Fermi liquid with spin-orbit coupling (SOC) is expected to support a new kind of collective 
modes: oscillations of magnetization in the absence of the magnetic field. We show that these modes 
are damped by the electron-electron interaction even in the limit of an infinitely long wavelength 
(5 = 0). The linewidth of the collective mode is on the order of /Ep, where A is a characteristic 
spin-orbit energy splitting and Ef is the Fermi energy. Such damping is in a stark contrast to 
known damping mechanisms of both charge and spin collective modes in the absence of SOC, all of 
which disappear at g = 0, and arises because none of the components of total spin is conserved in 
the presence of SOC. 


Electron systems with spin-orbit coupling (SOC) ex¬ 
hibit rich physics, some of which may have technological 
applicationSfii^ equally rich is the physics of cold-atom 
systems with synthetic SOC4ii^ Combining many-body 
interactions with broken SU(2) symmetry, one obtains 
a special-“chiral”-kind of Fermi liquid (FL)^“— that sup¬ 
ports a new type of collective modes, “chiral-spin waves”- 
oscillations of the spin density in zero magnetic field 

In the absence of SOC, electron-electron interaction 
(eei)^ does not affect certain properties of an electron 
system given that some symmetries are preserved. For 
example, the conductivity and cyclotron-resonance fre¬ 
quency of a Galilean-invariant system are not affected 
by eei; same is true for the de Haas-van Alphen (dHvA) 
frequency in an isotropic systent^^ and for the Farmer 
frequency in the presence of an SU(2)-symmetric in¬ 
teraction. Being a relativistic effect, SOC breaks both 
Galilean (but not necessarily rotational) invariance and 
SU(2) symmetry and thus lifts the protection ensured by 
these symmetries. As a result, several physical quanti¬ 
ties become dependent on eei. The list of such quantities 
includes the optical conductivity,— Drude weight,—, and 
frequencies of collective spin modes, which play the role 
of Farmer frequencies in zero magnetic field4i“— 

In this Fetter, we discuss another fundamentally new 
effect induced solely by SOC: intrinsic damping of collec¬ 
tive spin modes in the uniform (q = 0) limit. Interaction- 
induced damping of collective modes is not, by itself, 
a new effect. For example, plasmons in 3D)^ 2D,^i 
and ID (Ref. [H) electron systems, the Silin-Feggett 
(SF) mod o^^d^ in a partially spin-polarized FF )^^d^ and 
magnons in a ferromagnetic FL— are all damped by 
interaction processes involving excitations of multiple 
particle-hole pairs. (This mechanism is different from 
Fandau damping which involves only a single particle- 
hole pair: damping by multiple pairs occurs even outside 
the single-particle continuum.) However, Galilean invari¬ 
ance, in the case of charge modes;^ and conservation of 
the total spin component along the field ( 5 * 3 ), in the case 
of spin modeS)^ ensure that this kind of damping van¬ 
ishes at g = 0. We show here that this is not the case for 


electron systems with SOC. 
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FIG. 1: Left: Schematics of the Silin-Feggett mode in a par¬ 
tially spin-polarized FL. Right: The chiral-spin modes in a 
FL with Rashba spin-orbit coupling. The shaded regions de¬ 
note the particle-hole continua. As is the Larmor frequency. 
As is the quasiparticle Zeeman energy, and Amin/max is the 
lower/upper boundary of the continuum at g = 0. 
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FIG. 2: Top: The ladder (RPA) series for the spin suscepti¬ 
bility. The boxed wavy line is the static effective interaction 
Ux- Bottom: Diagrams contributing to damping of the collec¬ 
tive modes. The wavy line in diagrams a-e denotes a dynamic 
interaction, V'eft(D)- 

Two-dimensional (2D) electron systems with 
momentum-dependent SOC, e.g., of Rashba or Dres- 
selhaus types, bear certain similarity to a partially 
spin-polarized Fermi gas. The latter has a transverse SF 
mode in the spin sector (see Fig. [U left) )^^'^^ while the 
former has three (two transverse and one longitudinal) 
chiral-spin modes (Di... Da in Fig. [1] right)which 







2 


correspond to oscillations of the three components of 
magnetization. Indices 1 — 3 label the Cartesian system 
with the 3 axis along the normal to the plane of a 2D 
electron gas. (Although SL-mode has been studied 
previously in 3D, the same mode should occur in 2D 
as well.) Conservation of S 3 ensures that the frequency 
of the SL mode at <7 = 0 coincides with the Larmor 
frequency in the absence of eei. In the presence of SOC, 
none of the three spin components is conserved. As a 
result, one obtains three distinct modes with frequencies 
renormalized by eei. In addition, as we show here, 
these modes have finite linewidth which, in order of 
magnitude, is given by the inverse transport lifetime 
of a quasiparticle with energy equal to the spin-orbit 
splitting. That the SL mode at g = 0 is not affected by 
eei follows already from the exact equations of motion 
for magnetization.— Diagrammatically, this occurs due 
to a cancellation between the self-energy and vertex 
graphs for the spin susceptibility<^ Such a cancellation, 
however, does not occur in the presence of SOC. 

The single-particle Hamiltonian of a 2D system with 
SOC can be written as (we set h= 1) 

"Hfe = ao +Xa-f{k), (1) 

where /i is the chemical potential, CTq is the 2 x 2 unit ma¬ 
trix, <7 is the three-dimensional vector of Pauli matrices, 
A is the SOC constant, and f{k) = —f{—k) is a 2D vector 
that depends on the details of SOC; e.g., / = (fc 2 , —fci, 0) 
for linear Rashba SOC. The single-particle Green’s func¬ 
tion is given by 

GiK) ='^^sik)gsiK), Ds(fc) = i [cTo-b s%] , ( 2 ) 

S 

where gs{K) = {iko - ^ - sAj^/2 + g)-^, K = (ifcp, k), 
s = ±1 labels either spin projection or chirality, fj^ = 
CT-//I/I, and Aj: = 2A|/| is the spin-orbit splitting which, 
in general, depends not only on the magnitude but also 
on the direction of k. We will be primarily interested in 
the case of weak SOC, when Aj: <C Ep for any k. In 
what follows, we will be comparing the D 3 chiral-spin 
mode to the (2D) SL mode, as both modes are trans¬ 
verse to the SOC-induced/Zeeman magnetic field. The 
latter can be described by the same Hamiltonian with 
/ = (0,0, Ab/ 2A), where Ap = ggsB, g is the effective 
g—factor, /is is the Bohr magneton, and B is the mag¬ 
netic field chosen to be along the 3-direction. The orbital 
effect of the field is not considered here. 

Within the Random Phase Approximation (RPA), the 
spin susceptibility tensor is given by the ladder series in 
Fig. m where the boxed wavy line is a short-range inter¬ 
action, C/x, which mimics the exchange interaction in the 
spin channel. As shown in Ref. El the frequencies of the 
collective modes correspond to the roots of the equation 
Det(cro C) (To -b %-n°), where the elements of the 4x4 


spin-charge polarization matrix H® are given by 


n° (Q) = f Ty[aMK)<JjGiK + Q)], (3) 

JK 

where = TJ2koI ^ 0 . 1 . 2 ,3, and 0 corre¬ 

sponds to the charge component. At g = 0, all mixed 
spin-charge susceptibilities, Hq^ with j ^ 0 , vanish by 
charge conservation, while the matrix of spin suscepti¬ 
bilities can always be transformed to a diagonal form. 
In general, there are three spin modes, whose frequen¬ 
cies are found from the equations 1 -b = 0 , with 

i5i = (52 = 1 and i 53 = 1/2. The Green’s functions in 
n/j contain the self-energy parts. However, for the spe¬ 
cial case of Dx =constant, they drop out (see below) and, 
after analytic continuation (iqo —>■ D -b ii5), one obtains 


n°3(D) = 2i/ 


A2 

kp 

(n + iS)^ - Ai 

kp 


FS 


(4) 


where we have already assumed that SOC is weak in the 
sense specified above, v is the density of states per spin 
projection, (...) denotes averaging over the Fermi surface 
(FS), kp = kpk/k, and kp is the Fermi momentum in 
the absence of SOC. In general, A^ varies from Amin 
to Amax along the FS. The continuum of inter-subband 
particle-hole excitations, where Imn 33 7 b 0 , is confined 
to the interval Amin A O < Amax- At the boundaries 
of the continuum, ReHgg has square root singularities^ 
which guarantee a solution of the eigenmode equation 
1 -b t/xn 33/2 = 0 for D < Amin even at weak coupling. 
If SOC is isotropic, Amin = Amax = A, the two square- 
root singularities merge into a single pole at D = A, the 
inter-subband continuum shrinks to a single point, and 
the mode frequency is given by D 3 = A\/l — u, where 
u = 

Renormalization of the mode frequency by eei in the 
SOC case and the lack thereof in the SL case is an impor¬ 
tant difference, which we discuss now as it will help us 
to understand the differences in damping later on. This 
difference occurs because the Greens’ functions in the 
RPA series in Fig. [T] include chirality- or spin-dependent 
shifts in the chemical potential, which are given by the 
momentum- and frequency-independent parts of the self¬ 
energy. Each rung of the ladder diagram (HF) contains 
a difference of the self-energies 


(5Ss = Ss - =-Dx^ / {Bs,s' - B^s,s')gs'{P), 


(5) 


where Bg^s' is the matrix element for the transition 
s —>■ s'. Within a given rung, 5Es renormalizes the spin¬ 
splitting perturbation, be it SOC or the magnetic field. 
Since s = s' in the SL case, the self-energy of an electron 
with given spin is proportional to the number density of 
electrons with the same spin. Hence, JEg is proportional 
to magnetization [i5Ss = suAb /(1 —u)], and each rung of 
the diagram contains the renormalized Zeeman energy of 





3 


a quasiparticle, As = Ab/( 1 — u). The boundary of the 
continuum at g = 0 is shifted from Ab to A^ (cf. Fig.[Tl 
left) as can be seen, e.g., from the 11 component of the 

rung: 115^(11) = — . The Larmor theorem 

is effected via a cancellation between the self-energy and 
vertex contributions)^ when the rung is substituted into 
the eigenmode equation, the factor of 1 — rt cancels out 
and the frequency of the mode coincides with bare A b ■ 
The SOC case is different in that chirality, in contrast 
to spin, is not conserved by ee% and the sum over s' in 
Eq. ([5]) contains both the s = s' and s' = —s terms. 
For Ux = const, this implies that the self-energy of an 
electron with given chirality is proportional to the total 
number density, and thus AEg = Oj^ The vertex part 
is, however, non-zero. Therefore, there is no cancellation 
between the self-energy and vertex contributions, and the 
frequencies of the modes are renormalized by eei. If ^ 
const, one can sho w^^i^° that <5Es does not contain the 
zeroth angular harmonic of the interaction (which is why 
i5Ss = 0 for Ux = const), whereas the vertex contribution 
does. Thus, there is no cancellation between the two 
contributions in the general case as well. 

Damping of collective modes in the region of fre¬ 
quencies and momenta outside the single-particle contin¬ 
uum occurs via generation of multiple particle-hole pairs, 
which requires a dynamic interaction, e.^., a dynamically 
screened Coulomb potential. The self-energy of collec¬ 
tive modes is depicted diagrammatically in Fig. [2] a-e. 
The same set of diagrams has been encountered in the 
analysis of various two-particle correlation functions in 
the case of an RPA-type interaction4ir— Although the 
Aslamazov-Larkin (AL) diagrams d and e contain two 
wavy lines, they are of the same order in the bare cou¬ 
pling constant of the theory (the electron charge in our 
case), as diagrams a-c. However, the contribution of the 
AL diagrams to damping vanishes within the approxima¬ 
tions made in this work 4 Si 21 

In the case of the SL mode, renormalization of the 
transverse spin susceptibility (x_l ^ xii + X 22 ) by dia¬ 
grams a-c is given by 6 x±{Q) = -Hb + ^b+ ^c), 
where 


n„ = 


f g-{K) 9 +{K Q)Y._{K) -\- {a^^ —>• —A^^^ 

K 


Hf, 


He 


[ [gliK + Q)g+iK)^_iK + Q) 

Jk 


r r g.{K + Q)g+{K) 
Jk \ *9o + 


A 


kp 


-A 


kp 


)> 


[E+(iF)-E_(iF + Q)] 

( 6 ) 


± denote up/down spins, Q = {iqo,0), T,±{K) = 
— Jp g±{K P)Ves{P), and 14ff(P) is some dynamic in¬ 
teraction. In the last line of Eq. (0), we used the identity 


iqo ± Ab 

(7) 

and integrated over P. Because the denominator in 
Eq. d?]) does not depend on P, this last step produced the 
same self-energies, E±, as in diagrams a and b. Adding 
up the three lines of Eq. m , we arrive at 


Hq -|- Hf, -|- He 


2 Ab 
lo + 


(A+ — A_), 


( 8 ) 


where A± = fj^g^(K)E±(K). Recalling that t4ff(P) is 
real on the Matsubara axis and changing the variables as 
fco —>■ — fco and po —>■ —po, we find that A± = A^. Thus 
the frequency-independent prefactor in Eq. ([8]) , A+ — A_, 
is real. Continuing iqo to the real axis, we see that the 
imaginary part of 5xi. comes only from a resonance at the 
bare Larmor frequency, XI = Ab, which coincides with 
the pole of x_l in the RPA approximation. The only effect 
of the interaction processes represented by diagrams a- 
c is thus to renormalize the amplitude of the SL mode 
without either shifting its frequency or smearing it. 

For the case of SOC, it is convenient to consider 
renormalization of the out-of-plane spin susceptibility, 
<5X33(0) = (Tla + Hb -I- He), where now 


Ha 

Hb 

He 

Adrts 


i [gUK)g+{K + Q){E+_(iX) + S_+(iX)} + gl{K)g.{K + Q){^++(K) + E__(iX)}] , 

J 2 [d-i^ + Q)ff-i-(77){^H—(77 + Q) + E— ^.[K Q)} g\{K Q)g-{K){T,pp{K -I- Q) -I- E— [K Q)}] 


Af-+M+— J\f+-M+-+ Af-+M-++ \ 

} 

[9r{K + P) — gt{K P P + Q)] 1 + scos , Nrt = 9r{K)9t{K P Q)- 


\ iqo P iqo P “ ^k+p *90 - A 


(9) 
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Here, r, s = ± label the spin-split bands, (j)): de¬ 
pends on the azimuthal angle 9j: of k (and is equal 
to 9j: for linear Rashba SOC);^ and the partial 
self-energies are defined as T,rt{K) = — jpgr{K + 


P)Ves{P) I + t cos{4>i^ — : such that the total self¬ 
energies of the Rashba subbands are E+ = S++ -|- S_ 

and E_ = _-|- E_In the last line of Eq. m , we 

only used the identity ©• In contrast to the SL case, 
however, the spin-orbit splittings in the denominators of 


He depend on k+p, and thus integration over P does not, 
in general, produce the self-energies. This already tells 
us that, in general, the imaginary part of 6 x 33 cannot 


cancel out between the self-energy and vertex diagrams. 


However, there are two realistic approximations, 
namely, of a long-range interaction (p <C kp) and of weak 
SOC (A^ Ep), within which the momentum depen¬ 
dence of A^^- can be neglected. Assuming that these 
two conditions are satisfied, the vertex part can again be 
rewritten in terms of the partial self-energies. Even in 
this limit, however, there is no complete cancellation be¬ 
tween the self-energy and vertex diagrams. Namely, we 
find that Ha-t-Hfj-l-nc can be rewritten as Hu-t-HD, where 

nR = Ik {^+{K)gl{K) - E.{K)g^_iK)} and 


nD = - 


A 


kp 


K Qo 




{[E-{K + Q)- E+_(if)] g-{K)g+{K + Q) - (a,^^ ^ -A,-^) } 


( 10 ) 


The first term. Hr, has the same structure as in Eq. 
using the same arguments as before, we conclude that 
Hr does not contribute to damping. In contrast, the 
second term, Hd, does have, in general, an imaginary 
part at all frequencies, which means damping. To cal¬ 
culate Hd explicitly, one needs to specify the interac¬ 
tion, which we choose to be in the form of a dynami¬ 
cally screened Coulomb potential. Deferring the compu¬ 
tational details to Sec. IV of the Supplementary Material, 
we quote here only the final result for X 33 ; near the res¬ 
onance at H = H 3 , 


X3y (D) = {2vpl)-^A [Hi - (H + *r/2)2] 


A = 




( 11 ) 


where luq = r^Eplnr^ ^/IStt, = \/ 2 e^lvp is the cou¬ 
pling constant of the Coulomb interaction, and ^2 is a 

k p 

dimensionless form-factor which depends on the details 
of SOC; for isotropic SOC, 0 = 1. 

k p 

The damping rate T has an expected EL form. Notice 
though that the quasiparticle damping rate in 2D scales 
as Tqp oc In H, as opposed to just with a prefactor 
which does not depend on Being a gauge-invariant 
quantity, T contains the differences of the single-particle 
self-energies [see Eq. (HH)], while Tqp is related to the 
self-energy itself. The infrared singularity in the self¬ 
energy, which gives rise to the In H factor in Tqp, cancels 
out in r. As a result, T is on the order of the transport 
decay rate, which is much smaller than Tqp. The EL 
nature of the result for T indicates that it would not 
change substantially if, instead of a EL with Coulomb 
interaction, we would consider a EL of neutral particles 
with short-range interaction. The only change would be 


in loq which, for the case of a contact interaction with 
coupling [7, should be replaced by ~ {UvEpY. 

Since the frequencies of chiral-spin modes are propor¬ 
tional to the spin-orbit splitting, one might be tempted to 
conclude that it is better to look for these modes in mate¬ 
rials with strong SOC. Our result in Eq. (HU shows that 
the advantage of strong SOC has its limits. Indeed, the 
ratio of the linewidth to the mode frequency, 7 = E/A, 
scales as CA/Ep, where f and A are the appropriate 
angular averages of T and Aj^, correspondingly, and C 
a dimensionless prefactor. In a material with sufficiently 
strong eef one should expect that C ^ 1. If, in addi¬ 
tion, SOC is also strong (A ^ Ep), then 7^1 and the 
mode is overdamped. We emphasize that this effect is a 
unique feature of SOC; in contrast, the SL mode remains 
undamped (at q = 0) even if the Zeeman energy becomes 
comparable to the Fermi energy. 

In a particular model of the screened Coulomb poten¬ 
tial, the effect of damping appears to be rather weak. 
For isotropic SOC, Eq. (fTT|) yields C = r^lnr^ ^/127r. Us¬ 
ing parameters for an InGaAs/InAlAs quantum well, we 
then find 7 Ri 2 x 10 ~^ A/Ep for an electron number den¬ 
sity of 1.6 X 10^^ cm“^. On the other hand, chiral-spin 
waves are also damped by disorder via the Dyakonov- 
Perel’ mechanism;^ For a mobility of 2 x 10® cm^/V- s, 
damping due to disorder is stronger than that by eei by 
a factor of ten. One should not forget, however, that 
Eq. (fTT|) is valid only for Ts -C 1 and the actual numbers 
may differ from quoted above as increases. Although 
damping from disorder appears to be the dominant effect 
in solid-state systems, damping due to interaction should 
be dominant in (fermionic) cold-atom systems with syn¬ 
thetic SOC,-^ which have virtually no disorder. In this 
case, the interaction is short-ranged but, as we have al¬ 
ready mentioned, this should only affect the prefactor in 
Eq. HID. 

In conclusion, we showed that eei in the presence of 
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SOC not only gives rise to a new type of collective modes 
but also leads to their damping. This damping occurs 
even at 9 = 0 and its rate scales as the square of the spin- 
orbit splitting. This effect occurs because neither of the 
three components of magnetization is a good quantum 
number in the presence of SOC. 
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Supplementary Material 


I. THRESHOLD SINGULARITIES IN EQUATION (4) OF THE MAIN TEXT (MT) FOR GENERIC 

SPIN-ORBIT COUPLING (SOC) 


For a generic (but still weak compared to the Fermi energy SOC) the spin-orbit band splitting, is anisotropic. 

Suppose that varies from Amin to Amax along the Fermi surface (FS). Quite generally, the angular dependence 
of A^^^ near the extremal points is Amin ~ A_ -|- and Amax ~ A+ — (3^9“^. Near each of these extremal points, 
the diverging part of the angular integral in Eq. (4) of MT can be written as 


/ 

\ D2 _ A?, 

\ kp 



1 f de 

2 J 27: 

1 r dd 

2 J 27: 


A-i- 1 

(n-A++pif)^) 

A+ 1 


near maximum, 
near minimum. 


Note that the angular average is imaginary when A_ < D < A_|_, which correspond to the continuum of inter-subband 
particle-hole excitations. The singularity in the real part of the susceptibility at the lower end of the continuum ensures 
that the collective modes exist even for infinitesimally weak electron-electron interaction {eei). Eor an isotropic SOC, 
e.g., either for Rashba or Dresselhaus SOC, A+ = A_ = A and the two square-root branch cuts merge into a simple 
pole 1/(0 — A). 


H. STATIC PART OF THE SELF-ENERGY FOR MOMENTUM-DEPENDENT INTERACTION 


In the MT, we argued that the difference of the self-energies of chiral subbands does not contain the zeroth angular 
harmonic of the interaction potential. In this section, we prove this statement. The self-energy of a subband with 
chirality s = ±1 due to eei via a static but momentum-dependent potential, is equal to 

^s{K) = {5+(i^') +5-(^') + s[ff+(^') -5-(^')]cos(/gg,} (12) 


For the case of linear Rashba SOC, is the angle between k and k' (= for a general case, is some 

function of O^p such that / cos (f)^pd9^p = 0. We hrst consider the case of Rashba SOC in some detail and then 
point out the qualitative features that remain the same for arbitrary SOC. In addition, we also assume that Rashba 
SOC is weak (A <C vf), in which case the spin-orbit splitting can be approximated by its value projected on the FS 
in the absence of SOC: A « 2XkF- The difference of the self-energies of the chiral subbands, SUg = — E_s, is now 

given by 


S^s{K) = -s np (^sp + j^-nF (sp - 




where npie) is the Fermi function. Setting T = 0 and expanding the Fermi functions in A, we find 

dOz 


S'Ss{K) = sfA 


_ hhL rnsd— r/ - - 


(13) 


(14) 


where k'p = kpk'/k. To find the frequency of the collective mode, one needs to project k on the FS upon which SX3s 
becomes proportional to the first harmonic of the interaction potential. 

A similar argument can be made for an arbitrary (but still weak) SOC. The Green’s function in Eq. (2) of MT can 
be separated into symmetric (S) and asymmetric (A) parts: G = + G^. Any crystalline plane has at least a C 2 
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symmetry upon which k —>■ —k. This guarantees that A_j: = A^r and thus the subband Green’s functions, g±, are 

even on fc —>■ — fc. Then G® is even while is odd on this operation. Consider the self-energy in the spin basis (a 
2 x2 matrix) 


= -I^U^k-k'\GiK'), (15) 

which can be also decomposed into symmetric and asymmetric parts: E = E® -|- E^. The zeroth angular harmonic of 
the interaction potential, survives only in the part of the integral associated with G® and enters the diagonal 

elements of E as where n is the total number density. When E is transformed to the (diagonal) chiral basis, 

the terms remain on the diagonal and cancel out in the difference E+ — E_. 


III. ASLAMAZOV-LARKIN (AL) DIAGRAMS 

In the presence of a dynamic and long range interaction, the AL diagrams are of the same order as the self energy 
and vertex diagrams (See Refs. 29-32 of MT). In the absence of SOC and for a spin-invariant interaction, however, 
AL diagrams for the spin susceptibility vanish for a trivial reason: Tr[tTiGGG]= 0, {i = 1,2,3) because G oc ctq. 
In the presence of the magnetic field (along the 3-axis), G = aag + ha^. The relevant quantity in this case in the 
transverse spin susceptibility, which contains the trace Tr[criGGG] with i = 1,2. Again, this trace is equal to zero. 
The vanishing of the AL diagrams in both these cases is a consequence of conservation of either total spin (in the 
former case) or its component along the magnetic field (in the latter case). In the presence of SOC, this reasoning 
does not apply because G now has off-diagonal elements. However, some general statements about the AL diagrams 
can still be made for q = 0 case. 



FIG. 3: Detailed AL diagrams from Fig. 2 of MT 


For convenience, we present the AL diagrams from Fig. 2, d and e of MT here in Fig. [3] It suffices to consider only 
the diagonal components of the spin susceptibility, Xmj with f = 1... 3, in which case the diagrams read 

n-. = f f Tr[a,G{K)G{K + P)G(K + Q)]VpVp.Q f Tr [a,G(Ai + Q)G(iLi + P)G(Ai)], 

JPJK JKi 

n-. = f f TT[a,GiK + Q)GiK + Q-P)GiK)]VpVp.Q f Tr [u,G(iLi + Q)G(Ai + P)G(Ai)]. 

JpJK JKi 

(16) 

Here, Q has only the temporal component: Q = {iqo, (1). We can define 


hi{P,Q) = [ Tr [(TjG_r-+qGk+pG_r-] , 

Jk 

hi{P,Q)= [ Tr [cTjGifGif+pGic+Q], 

JK 

fi{P,Q)= [ Pt^[o'iGk+qGk+q-pGk] ■ 

Jk 

Using the matrix structure of G{K), we observe that 

hi,2{P,Q) = +hi,2{P.Qy. h^iP, Q) = -hiP^Q) a.iid /,(P, Q) =Q) for z = l,2,3. 


(17) 


(18) 
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With the help of these properties, the AL diagrams can now be written as 


nf, = / h^{P,Q)h^{P,Q)VpVp-Q = / h?(P,Q)ypVp_Q fori = 1,2 and 


= - J h‘^iP,Q)VpVp_Q for i = 3, 


n?. = f{P, Q)h^iP, Q)VpVp-q = -J^ |h.(P, Q)\^VpVq-p, 


and thus their sum is reduced to 


nf. + nj. = {hj{P, Q) - MP, Q)p) VpVp.Q, for * = l, 2 and 
= - J^{h^,iP,Q) + MP,Q)f)VpVp-Q, for i = 3. 


(19) 

( 20 ) 


( 21 ) 


We now study two particular cases-that of Rashba SOC and of combined Rashba-Dresselhaus SOC- and show that 
the AL diagrams do not contribute to damping in either of these two cases. 

Rashba SOC: For i = 3, 


/i3(-P,Q) = -^ X! b(a-c)J sm{0i:-0i:_^^)ga{K+ Q)gb{K + P)gc{K). 

a,b,c=± ^ 


( 22 ) 


Since Q has no spatial component, we can choose the direction of p as a reference. On reflecting vector k about this 
direction, both 0^ and 0^_^_p change signs and so does the factor of sin(...). On the other hand, since the energy 
spectrum is isotropic, each of the three subband Green’s function in the formula above is a function of the magnitude 
of the corresponding momentum. It is then easy to see that all the three Green’s function are even on reflection 
ki ^ ki, k 2 ^ —k 2 - Therefore, the integral over 0^ vanishes, and the AL diagrams give no contribution to X 33 . 

The same reasoning applies for i = 1. In this case. 


hi(P,Q) = - ^ 

a,fe,c—± 


K >- 


(a + c) sindg + a&csin(26»g + %^-) + 6 sin%^- ga{K + Q)gb{K + P)gc{K) (23) 


contains again an angular factor, which is odd on reflection, and a combination of the Green’s functions, which is 
even on reflection. Therefore, hi{P, Q) = 0 and there is no contribution to xii from the AL diagrams either. 

The same argument cannot be used for * = 2, because the angular factor in 


h 2 {P,Q) = -^ ^ J {a+ c) cos 0j: + abccos ( 20 + 0+ boos0i^_^_p) ga{K + Q)gb{K + P)gc{K). (24) 

a,b,c—± ^ 


is even on reflection. Nevertheless, in-plane rotational symmetry of the Rashba Hamiltonian ensures that xii = X 22 - 
Therefore, one does not need to consider X 22 separately: if there is no damping of resonance associated with xii-, the 
same is true for X 22 - 

Combined Rashba and Dresselhaus SOC: In the presence of both Rashba and Dresselhaus SOG, the only change in 
Eqs. ( (|22II24I) 1 is that 0j: needs to replaced by cj)^, as defined after Eq. (fl^ . It is convenient to choose the Gartesian 
system rotated by 7 r /4 compared to the conventional one. In such a system, the combined Rashba plus Dresselhaus 
Hamiltonian is described by in Eq. (1) of MT with A/(fc) = ([/3 -|- q;]A: 2 , [fd — a]/ci, 0) and 


cos^j = 




(a - 13) cos 

(/3 -b a) sin 0^ 

% 


(25) 


where Ag = -b — 2a/3 cos 20^. Note that = tan ^ ^^^tan^g^ The spin-orbit splitting = 2kpA^ 
is even on reflection about the fci-axis. As discussed in MT one can approximate A^_|_^ by for a long-range 
interaction; therefore, the products of the Green’s functions are still even on reflection. Because of the relation 
between and 0^, the same symmetries that led to the vanishing of hi and h^ in the case with only Rashba SOC 
ensure that hi and h^ vanish is the case when both Rashba and Dresselhaus SOC are present as well. 










IV. DERIVATION OF EQUATION (11) OF THE MAIN TEXT 


In this section, we present the evaluation of IId from Eq. (10) of the MT. For convenience, we present this equation 
below in an expanded form: 


Hd 




+ Q) - S+_(A)] g-{K)g+{K + Q) - [S+_(K + Q) 


E__(A)]g+(if)g_(A + Q)}, 


(26) 


where Y,ab{K) = — Jpga{K + P)Ves{P) 1 + 6cos(0^ — 0^^^) and the angle is defined after Eq. 

is already proportional to we can neglect the effect of SOC on the interaction potential and 
dynamically screened Coulomb potential in 2D: 


(13. Since IId 
take it to be a 


V,s{P) 


P + K 


27re^ 


(1 , IPol ^ ’ 

\ yJvl+PpP"^ J 


(27) 


where k = \/2TskF is the inverse Thomas-Fermi screening length and = v^e^/uF- The momentum transfer p 
is decomposed into two components: along the normal to the FS at point k (p||) and along the tangent (p_l). We 
assume and then verify that the quasiparticle interaction is determined by processes with typical p|| on the order 
of po/vf, whereas typical p± are on the order of k. In turn, typical po are on the order of the external frequency 
which, in our case, is fixed by spin-orbit splitting, A, averaged over the FS. We choose to work in a realistic regime 
of A <C vfk <C Ef- [Although the form of the potential in Eq. (I27p is, strictly speaking, valid for rg <C 1 which 
implies that k <C /cf, these strong inequalities are never satisfied in real materials. Therefore, the interval of energies 
in between vfk and Ef is never wide enough to consider a possibility of A being within this interval.] Since 


bill ~ bol/^^F < bu| (28) 

in this regime, we can simplify the Landau-damping term in Eq. (1271) to bol/nFbj-l and also expand the denominator of 
this equation in this parameter. Subtracting off the static screened potential which gives no contribution to damping, 
we obtain for the dynamical (po-dependent) part of the interaction 


rrdyn 


(P) = 


27re^ 


bo I 


\p±\+kJ vf \p±\' 


(29) 


Thanks to the inequality (E51) . the integrals over the fermionic (K) and bosonic (P) momenta become separable. Let 
us pick one of the terms in Eq. (l26l) : 



+ Q)g_{K)g+{K + Q) 


[ [ U^rg.iK + Q + P)g_{K)g+{K + Q) 
Jk jp 


1 — cos 




■ (30) 


Since it is expected that typical |p_l| ^ k kF while k is expected to be near the FS, the angular factor can be 

2 

expanded as I — cos{cj)^ — (/'j;_|_^) ~ where is a form-factor that depends on the details of the SOC. 

For example, in the case of linear Rashba SOC, (pj: coincides with the angle 9^ that k makes with, e.g., the I-axis. 

Also, due to rotational symmetry, the result of the integral over P does not depend on the direction of k and thus we 

2 

can choose = 0 in 1 — cos(9^ — 9^^^). An expansion of 1 — cos(9^p^) then yields which means that = I for 

a system with linear Rashba SOC. Notice that the condition |p||| <C bj-l implies that p is almost perpendicular to k, 
i.e., that bp — ~ tt/2 . 

In the presence of both Rashba and Dresselhaus SOC, pj: is given by Eq. dUD. Expanding in p and using that 
bp - « 7r/2, we get 


p 2 

1 - cos(bp - bp+ff) ~ 7 ^ 


1 -b 


7 


1 — 7 cos 207 


V. 2a/3 

where 7 = ^5 - ^ < 1 

+ j3^ 


(31) 


Thus ^1 = 1-1- system. We avoid the special case of |a| = |/3|, when 7 = 1 and has a pole: it 

can be shown that collective modes are completely covered by the continuum in this case.'^ 
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Next, we write g-{K)g+{K + Q) as [g-{K) — g+{K + Q)]/{iqo — Ag) and integrate the products of two Green’s 
functions over fcg and p||; for one of the products, we obtain 




Wo 


StT ivn — ‘ 


- 


bo I 
2 vp 


ipo - vfP\\ — 

The rest of the integrals are evaluated in the similar manner. Collecting all terms, we arrive at 
// Ar fr 

nD = - 




2 v 


rkp 7 poo 7 

dp± 2 / ^PO 


9o 




Vpk^ 


F JO 


-j^Pl -^UTiP) {\qo+Po\ + bo - Pol - 2|po|} . 


(32) 


(33) 


The factor in {... } constraints the range of integration over po to (0, qq). This leads to 


n 


D 



pkp 


dpo{qo-Po)Po I dp± 


P± 


{p± + nf 


(34) 


To logarithmic accuracy, the integral over pj_ gives Inkp / — Inrg while the integral over po gives Qq/G. Note that 

typical Po ^ qo while typical p_L n (in the logarithmic sense), and thus our initial assumption is verified. 

After analytic continuation, we get 


n^) = —iv 




02 - A2 

'v kp 



(35) 


where iJq = rgAplnr^ ^/127r and = \f2^jvp- This is the correction to Ilgg that is responsible for damping. 
Recalling that the physical susceptibility is given by X 33 = + ^ 1133 ) 1133 differs from Ilgg due to the 

damping correction. To account for this correction, we may write 


1133 = 2 v 


Al 

k p ^k p 


A 2^-(0 + * r / 2 )^ 


FS 


(36) 


where T = 4^ ' Using the fact that O 3 is the pole in X 33 when T = 0, we can eliminate the [4 dependence. 

Writing A| — (O + iT/2)^ as A| — 0§ + 0§ — (O + iT/2)^ and expanding in 0§ — (O + iT/2)^ near the resonance, 
we obtain Eq. (11) of MT. 
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